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The nature of the ASCA/ INTEGRAL source 
AX J183039-1002: a new Compton-thick AGN? 
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ABSTRACT 

We report on the identification of the X/soft 7-ray source AX J183039-1002 detected 
with ASCA and INTEGRAL/IBIS. The source, which has an observed 20-100 keV 
flux of ~ 8.6 x 10 -11 erg cm s _1 , is inside a diffuse radio supernova remnant (SNR) 
and is spatially coincident with a compact radio source. We analyzed archival Chandra 
and XMM-Newton observations in order to identify the ASC A/INTEGRAL source. 
A point-like Chandra X-ray object was found to be positionally coincident with the 
compact radio source and within the error circle of the ASCA and INTEGRAL sources. 
Although the association of a compact radio/X-ray source with a radio supernova 
remnant could be indicative of a pulsar wind nebula (PWN), the XMM-Newton X- 
ray spectrum is compatible with an absorbed, Seyfert-2 like AGN, since it provides 
evidence for an iron emission line of ~ 1 keV equivalent width; furthermore the X-ray 
source spectrum is similar to that of other Compton thick AGN where the <2 keV 
data are associated to a warm reflector and the >10 keV one to a cold reflector. 
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1 INTRODUCTION 

A key strategic objective of the INTEGRAL [International 
Gamma Ray Astrophysics Laboratory) mission is a high en- 
ergy survey of the Galactic plane and Centre (Winkler et al. 
2003) in primis and of all the sky as a by-product of pointing 
observations. This makes use of the unique imaging capabil- 
ity of IBIS {Imager on Board Integral Satellite; Ubertini et 
al. 2003) , which allows the detection of sources at the mCrab 
level, with an angular resolution of 12' and a point source 
location accuracy (PSLA) of typically l'-3' within a large 
(29° x 29°) field of view. So far, several surveys produced 
from the IBIS/ISGRI data have been reported, the most 
complete being those of Bird et al. (2007) and Krivonos et 
al. (2007); both reported around 400 sources down to a flux 
level of a few mCrab above ~20 keV. A significant fraction 
(25-30%) of the objects in these surveys have no obvious 
counterpart at other wavelengths and therefore cannot yet 
be associated with any known class of high energy emitting 
objects. Searching for counterparts of these new sources is of 
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course a primary objective of the survey work but it is made 
very difficult by the relatively large INTEGRAL error boxes. 
Cross correlations with catalogues in other wavebands can 
be used as a useful tool with which to restrict the positional 
uncertainty of the objects detected by IBIS and therefore to 
facilitate the identification process. Observations at softer 
X-ray energies, where the positional accuracy is much bet- 
ter, has already proved an invaluable aid in the identification 
and classification process (Stephen et al. 2006; Masetti et al. 
2008 and references therein), but the use of data at other 
wavebands, for example the radio, is also useful to help iden- 
tify peculiar and interesting objects. 

In this Letter we report on the unusual nature of 
the ASCA/ INTEGRAL source AX J183039-1002, located 
within diffuse radio emission and likely associated to a com- 
pact radio source. We use unpublished Chandra and XMM- 
Newton observations to localize and identify the X-ray coun- 
terpart of this high energy emitter and to characterize its 
broad-band spectrum. We provide arguments for the identi- 
fication of this X/gamma-ray source with either a new com- 
posite Supernova/PWN system or a background AGN; this 
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second option is however the more convincing on the basis 
of the source spectral characteristics. 



2 THE ASCA/INTEGRAL SOURCE 
AX J183039-1002 

AX J183039-1002 was first reported by Sugizaki et al. 
(200f) as an unidentified and faint X-ray source detected 
by ASCA during the Galactic Plane Survey, ft was located 
at RA(J2000) = f8 h 30 m 39 s and Dec(J2000) = -f0°02'42" 
with a positional uncertainty of ~3'. Its soft X-ray spectrum 
is well described (\ 2 /v — 4.9/10) by an absorbed flat power 
law with T = 0.04+° gg, iV H = (3.1^) x fO 22 cm" 2 and an 
unabsorbed 0.7-10 keV flux of ~2.3 x 10" 12 erg cm" 2 s" 1 . 

AX J183039-1002 was then reported as an IBIS source 
both by Bird et al. (2007) and Krivonos et al. (2007); here 
we use data collected in the 3rd IBIS survey (Bird et al. 
2007), which consists of all exposures from the beginning of 
the mission (November 2002) up to April 2006. The total 
exposure on this region is ~1.78 Ms. ISGRI images for each 
available pointing were generated in various energy bands 
using the ISDC offline scientific analysis software version 5.1 
(OSA 5.1; Goldwurm et al. 2003). The individual images 
were then combined to produce mosaics of the sky region 
of interest here to enhance the detection significance using 
the system described in detail by Bird et al. (2004, 2007). A 
clear excess is observed in the IBIS map with a significance of 
~8cr at a position corresponding to RA(J2000) = 18 h 30 m 39 ! :8 
and Dec(J2000)=-10°00'25'.'2 and with an associated uncer- 
tainty of 4.8' (90% confidence level). Our position is com- 
patible with that reported by Krivonos et al. (2007) (note 
that their 2.1' error radius corresponds to 68% confidence 
level). Figure 1 shows the 20-40 keV image of the region 
surrounding AX J183039-1002 with superimposed both the 
INTEGRAL and ASCA error circles. The IBIS/ISGRI spec- 
trum was obtained following usual procedures: fluxes were 
extracted from the location of the source in 10 narrow energy 
bands over the 17-100 keV range for all available pointings; 
a spectral pha file was then made by taking the weighted 
mean of the light curve obtained in each band. An appropri- 
ately re-binned rmf file was also produced from the standard 
IBIS spectral response file to match the pha file energy bins. 
Here and in the following, spectral analysis was performed 
with XSPEC v.11.3.2 package and errors are quoted at 90% 
confidence level for one interesting parameter (A^ 2 = 2.71). 
A simple power law provides a good fit to the IBIS data 
(X 2 /^ = 7.6/8) and a photon index V = 3.llg;? (much 
steeper than the ASCA one) combined to an observed 20- 
100 keV flux of 8.6 x 10" 11 erg cm" 2 s" 1 . 



3 RADIO COUNTERPART 

The Multi-Array Galactic Plane Imaging Survey (MAG- 
PIS, Helfand et al., 2006) has mapped the region of in- 
terest here with a 20 cm sensitivity of 1-2 mjy and de- 
tected a compact source surrounded by diffuse extended 
emission (see Figure 2 which is a MAGPIS extract). The 
compact source (G21. 63201-0. 00682) which is located at 
RA(J2000) = 18 h 30 m 3S!29 and Dec(J2000)=-10°02'46'.'4 (5" 
uncertainty) has an integrated 20 cm flux density of 5.13 ± 




Figure 1. IBIS/ISGRI 20-40 keV image of AX J183039-1002. 
The green and white circles correspond to the INTEGRAL/YBIS 
and ASCA/GIS error boxes. 



0.26 mjy, and is also detected at 6 cm with an integrated 
flux density of 1.89 ± 0.18 mjy (White, Becker & Helfand, 
2005). The diffuse emission (G21. 6417+0.0000) around this 
discrete source has an integrated flux of 1.54 Jy at 20 cm 
and a diameter of 2.8'. 

Comparison with mid-infrared images obtained with 
MSX (Price et al., 2001) allows us to discriminate between 
thermal and non-thermal sources: indeed high values of the 
radio-to-IR flux ratio (20 cm versus 20 /im, where the source 
is generally not detected in the MSX) generally indicate non 
thermal radio emission such as is produced in a supernova 
remnant (SNR, see, e.g., Fiirst et al. 1987, Green 1989). This 
evidence together with the morphology of the diffuse radio 
emission (in this case almost a complete shell structure) and 
the detection of the source in the 90 cm images, led Helfand 
and co-workers to propose G21. 6417+0. 0000 as a high prob- 
ability SNR candidate. 

The association of a hard X-ray source coincident with 
a radio point source and surrounded by a very low surface 
brightness radio shell is suggestive of a situation where a 
SNR harbors a classic young radio pulsar possibly together 
with its wind nebula (PWN). A displacement of the compact 
source with respect to the centre of the shell structure is not 
unusual, since many pulsars are born with high kick veloc- 
ities which allow them to propagate through the shocked 
ejecta in the SNR interior (Gaensler and Slane 2006). Al- 
though the distance and age of G21. 6417+0. 000 are not well 
known, it is clear from its small angular size that it's rela- 
tively small and young. If the compact source is a PWN, 
therefore, its offset from the centre of the shell would re- 
quire an unusually large kick velocity. A way to know if our 
discrete radio source is a PWN is by means of its radio spec- 
trum, which is expected to have a flat energy index 0.3 ^ 
a ^ 0.1 (defined by S v oc v~ a , Chevalier 1998): comparison 
between 20 and 6 cm data provides an a value of 0.9, con- 
siderably steeper than what generally found in PWN. Also, 
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Figure 2. MAGPIS (20 cm) image of AX J183039-1002, showing 
both the diffuse emission and the compact radio source. 



the ratio between radio and X-ray core luminosities is higher 
(~2.3 x 10 4 ) than usually found in PWNe (Becker & Helfand, 
1987), indicating either a highly unusual PWN or an alter- 
native nature for AX J183039-1002, for example that of a 
background AGN. Indeed the radio spectrum of 0.9 is typical 
of an active galaxy, although this scenario is disregarded at 
first glance for the low probability of a chance alignment be- 
tween an AGN and the SNR. However a more in depth anal- 
ysis indicates that the probability of such an association is 
small but not insignificant. The unabsorbed 2-10 keV X-ray 
flux of the likely counterpart of AX J183039-1002 is about 
2xl0~ 12 erg cm" 2 s -1 , and from the LogA r -LogS' relation 
(La Franca et al., 2005) we estimate a density for such bright 
AGNs of about 0.1 sources per square degree. The SNR 
has a radius of about 2.8', i.e. an area of 6.15 arcmin 2 . We 
have then a mean value of about 2xl0~ 4 AGNs inside the 
SNR. However, we have observed in all the sky about 200- 
300 SNRs (for example Green's catalog reports 265 sources, 
see http://www.mrao.cam.ac.uk/surveys/snrs/). Assum- 
ing that all of these have the same area of G21. 6517+0. 000, 
we can next estimate the probability of finding at least 
one chance association with an AGN to be about 5%. This 
probability is small but not insignificant. Moreover, many 
SNRs have a radius much larger than 2.8', therefore the 
"true" chance association probability should be considerably 
greater than ~5%. 



4 CHANDRA OBSERVATION 

On February 19, 2006 Chandra carried out an observation of 
the sky region containing AX J183039-1002. Data reduction 
was performed with CIAO v. 3. 4 and CALDB v. 3. 2. 4 to ap- 
ply the latest gain corrections. Subsequent filtering on event 
grade and exclusion of high background times resulted in a 
total exposure of 19.2 ks. In order to locate the source of 
hard X/gamma-ray emission associated to the ASCA/IBIS 
object, we concentrate the Chandra image analysis to the 3- 
10 keV band: within the smaller ASCA error box we detect 



Table 1. 


Chandra X-ray 


detections. 






Source 


RA (J2000) 1 " 


Dec (J2000)t 


S/N 


Count Ratet 


Nl 


18 30 38.3 


-10 02 47.1 


27 


38.90 ± 1.50 


N2 


18 30 37.1 


-10 02 45.4 


4.7 


1.25 ±0.28 



t la uncertainty on Chandra positions is 0.6". 

t In units of 10 — 3 counts s — 1 in the 3-10 keV energy band. 



2 objects above a 3cr detection threshold; their positions, de- 
tection significance and count rates are reported in Table 1. 

Figure 3 is a cut-out of the 3-10 keV Chandra image, 
showing the 2 objects detected within the ASCA error box, 
together with the location of the compact MAGPIS radio 
source G21.63201-0. 00682. It is evident from the figure that 
the Chandra observation further confirms the association 
of source Nl with the compact MAGPIS radio object and, 
purely for its brightness, also with the ASCA/IBIS source. 
The reduced positional uncertainty does not contain any 
optical/infrared counterpart for this Chandra source in the 
HEASARC archive and even in the DSS infrared image (see 
Sect. 6 for a discussion on the expected optical extinction). 

To look for the possible presence of a PWN around 
source Nl, we have also compared the Chandra PSF against 
that of the source but have found no evidence for it: object 
Nl is seen as point-like by Chandra implying that a PWN 
is not present or, if present, is extremely compact. 

The CIAO script psextract was then used to generate 
the spectrum, with appropriate background and response 
files, of both source Nl and N2; spectra were extracted using 
a radius of 5'.' 6, while background files were generated using 
a circular region of radius 100". With a count rate of 0.119 
counts/frame, the pile-up fraction is insignificant at less than 
2%. In all our fitting procedures we used a Galactic column 
density, which in the direction of AX J183039-1002 is 1.52 
x 10 22 cm" 2 (Dickey & Lockman, 1990). 

The spectrum of source N2 is very soft with very few 
counts above 3 keV; when compared to the IBIS spectrum 
we have to assume a cross-calibration constant ^ 4x 10 4 , im- 
plying that the extrapolation of the Chandra data predicts 
an IBIS flux far below that observed, a further indication 
that this source is not the counterpart of the ASCA/IBIS 
object. 

The spectrum of source Nl is instead very hard as a sim- 
ple power law model provides a flat spectrum (P = —0.86) 
and an unacceptable fit (x 2 / u = 87.8/45). Residuals to this 
model show some curvature possibly due to intrinsic ab- 
sorption. Addition of this extra component provides a fit 
improvement {\ 2 /v — 40/45), which is significant at the 
99.99% confidence level, a steeper spectrum (r = 1.01 jig 38 )> 
a column density of A H = (ll^l 3 ,;®) x 10 22 cm -2 , and an 
unabsorbed 0.7-10 keV flux of (3.2l ( > 4 ) x 10" 12 erg cm -2 
s _1 . The Chandra spectrum is compatible within errors with 
the ASCA one. The spectrum is quite flat (as found previ- 
ously) and the column density exceeds largely the Galactic 
A^h value: it points to a source able to emit from a few keV 
up to ~100 keV, constant in time and, as underlined before, 
able to produce emission also at radio frequencies. 
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AX J 183030-1 002 - XMM-Newton absorbed powerlaw 
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Figure 3. Cut-out of the Chandra 3—10 keV band image showing 

the 2 objects detected within the ASCA error box, together with Figure 4. XMM-Newton spectrum of AX J183039-1002 fitted 

the XMM-Newton 1.06" lcr uncertainty as well as the location of with an absorbed power law: note the systematic residuals, 

the compact MAGPIS radio object. 



5 XMM-NEWTON OBSERVATION 

The Chandra Nl source has also been detected during the 
XMM Galactic Plane survey (Hands et al., 2004) as XGPS- 
I J183038-100248. It is also reported in the 2nd Serendip- 
itous Source Catalogue as 2XMM J183038. 2-100246 (Wat- 
son et al., 2007, http://xmmssc-www.star.le.ac.uk); the 
source 0.2-12 keV flux is ~2xl0~ 12 erg cm -2 s _1 (i.e. 
compatible with the ASCA /Chandra flux), with almost all 
counts coming from the hardest X-ray band (4.5-12 keV), 
again an indication of a flat 2-10 keV spectrum. 

We have then analyzed the spectrum of 
2XMM J183038. 2-100246, using its pipelined and re- 
duced 2XMM catalogue data products available at the 
HEADAS archivffl This observation was performed on 
September 13, 2003, for a total of 4335 seconds of exposure 
time with the EPIC-pn instrument. 

The source was constant during the entire observation 
and its spectrum is poorly described by a simple absorbed 
power law {x 2 /v = 62/19), at variance with that obtained by 
ASCA and Chandra. Excess counts are visible at low energies 
and at around 6.4 keV (Figure|4|. A good description of the 
data (x 2 / l/ = 12/14) is achieved with a model consisting 
of a blackbody (kT = 0.20 ± 0.02 keV), an absorbed power 
law (JV H = (8.8 ± 4.0) x 10 22 cm -2 ; F = 0.6 ± 0.5) and an 
iron K a line. Both these two extra features are required by 
XMM-Newton data, with a significance of 99.99% for the 
blackbody component and 99.6% for the line. This line is 
very little redshifted (z ^ 0.08) and has an equivalent width 
of ~1 keV (a — 0.4 ± 0.2 keV). Considering the errors on 
the line width and normalization given by the low statistics 
of the data, the EW could be significantly lower than 1 keV. 
If present, the broad line may be due to a blend of narrow 
lines (as often seen in Compton-thick AGNs, e.g. Matt et al. 
2004). Note that Chandra does not clearly detect the line. 
This is not surprising, since the effective area of ACIS in the 
Fe-K band is so small. 

Finally, we combined the XMM-Newton data with the 

http : //heasar c . gsf c . nasa . gov/W3Browse/xmm-newton/ 
xmmssc .html 



INTEGRAL one to obtain a broad-band spectrum over the 
1-200 keV energy band. The model used to fit the XMM- 
Newton spectrum gives residuals in the IBIS data suggest- 
ing the presence of some curvature or a cut-off above 10 
keV, which is immediately obvious also from the compar- 
ison between the X and 7-ray photon indices. To describe 
the broad-band source spectrum and in view of the low sta- 
tistical quality of the data we have therefore replaced the 
absorbed power law component with other simple models, 
like an absorbed broken power law or an absorbed cut-off 
power law; we have further fixed the cross-calibration con- 
stant C between XMM-Newton and IBIS to 1 in view of the 
non-variability of the source flux in the keV band given by 
comparing the ASCA, Chandra and XMM-Newton fluxes. 
We emphasize that the fit does not change considerably leav- 
ing C free to vary and in any case the large uncertainty on 
this parameter reaches as a lower boundary (7 = 1. The first 
model gives the same parameters found when using XMM- 
Newton data alone, a break at 26 ± 5 keV where the hard 
power law steepens to V > 2.8 (x 2 /v = 12/21). This fit is 
shown in Figure [5] The second model, apart from the other 
parameters that don't change significantly, has F = 0.0 ±0.5 
with a cut-off energy at 14 ± 5 keV (x*/v = 18/22). 

The two fits are equally acceptable and can be used 
to describe in a phenomenological way the source spec- 
trum, i.e a low-energy blackbody together with an absorbed 
flat power law steepening at around 10-30 keV and a low- 
redshift iron line. 



6 WHAT TYPE OF SOURCE IS 
AX J183039-1002? 

AX J183039-1002 is certainly a broad-band emitter over the 
1-100 keV band; it is also a weak and compact radio source 
sitting inside or behind a possible supernova remnant. No 
optical source is present within the restricted Chandra error 
box down to a limiting B magnitude of ~20-21; this suggests 
either that the source is weak or heavily obscured. Outside 
the Chandra positional uncertainty we find a DENIS and 
a 2MASS object at 1.1" and 2.6" angular distance, respec- 
tively. Both have similar near infrared magnitudes (J = 15.8 
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AX J1 83039-1 002 - XMM-Newton and INTEGRAL-IBIS/ISGRI 
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Figure 5. Broad-band spectrum of AX J183039-1002 fitted with 
a low energy blackbody, a redshifted iron line and an absorbed 
broken power law: XMM-Newton data to the left and IBIS data 
to the right. 



and K = 13) and compatible RA and Dec within the respec- 
tive positional errors so that they can be assumed to be the 
same source. While the DENIS position is compatible within 
errors with the Chandra one, the 2MASS location is too off- 
set from the X-ray one to claim a match; the positional 
accuracy of the 2MASS catalogue is ^1", which combined 
with the Chandra uncertainties provides at most an angular 
distance of 1.6", smaller than the observed one. The mis- 
match is however not too big and could be explained by a 
difference in object centroids, due for example to the source 
shape (as expected in an extended source like a galaxy) or 
to the source proper motion (as typical of a galactic source). 
We therefore assume that the DENIS/2MASS source could 
be the near-infrared counterpart of AX J183039-1002. 

The overall morphology of the source suggests two pos- 
sible scenarios: a PWN associated to the SNR or a back- 
ground AGN. In the first case the source would be a very 
unusual pulsar/PWN system from the radio and X/soft 
gamma-ray spectral characteristics. For example, the soft 
X-ray spectral index of a PWN is usually around 1.5-2 (Kar- 
galtsev & Pavlov 2008), softer than the ASCA and Chandra 
spectra. Furthermore, the photon index F > 2.8 observed 
in the IBIS band is significantly higher than the usual val- 
ues for the INTEGRAL PWNe, that cluster around F ~ 2.2 
(Dean et al. 2008). The alternative possibility that we are 
dealing with a magnetar is disregarded on the basis of the 
stable 2-10 keV flux over the years and the shape of the 
high energy emission (typically IBIS sees a hard tail in this 
type of sources, Kuiper et al. 2006). 

In the second, more likely, case we would deal with 
a heavily absorbed AGN (a type-2 Seyfert?) with an ex- 
tremely low exponential cut-off energy. The radio flux and 
spectral index would be compatible with the range of val- 
ues found in other INTEGRAL selected AGN, which are 
almost invariably associated with radio sources (Mushotzky 
2004). In this case, the spatial coincidence between a back- 
ground AGN and a SNR is fortuitous. But the more con- 
clusive evidence for the association with a heavily absorbed 
AGN comes from the XMM-Newton spectrum, that shows a 
strong iron line and is reminescent of a Compton-thick AGN. 



We have phenomenologically fitted the source spectrum with 
simple models, like a blackbody and broken power law. From 
a comparison with other similar sources in the literature (e.g. 
NGC 1068, Matt et al. 1997; see also Maiolino et al. 1998) we 
can associate the low-energy blackbody-like component with 
a warm reflector, and the high-energy component as a cold 
reflector coming from the inner torus. The ultimate solution 
to this puzzle will only come from detailed radio observa- 
tions: hopefully it will not be long before AX J183039-1002 
is definitely properly classified. If the heavily absorbed na- 
ture of this source is confirmed, this brings to 5 the number 
of candidate Compton-thick AGN in the third IBIS cata- 
logue, besides the 7 confirmed sources, i.e. about 20% of the 
total sample of type-2 active galaxies so far detected (Bird 
et al. 2007). 
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